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ROAD FRICTION TESTS IN ISRAEL 


M. PELEG 
Technion—Israel Institute of Technology, Haifa 


ABSTRACT 


Twenty-six points on the main roads of Israel were examined with regard to their friction 
properties. Using the Leroux “‘rugosimeter’’, the percentage was determined of kinetic energy 
of a pendulum absorbed by friction between its rubber lining and the road surface under 
dry and wet conditions. The safety coefficient of the road was calculated and the category 
of allowable speed determined. 

Observations have shown that for most of the roads tested only low and moderate 
speed driving is safe. It was found that the maximum aggregate size to be used for asphalt 
wearing courses, in order to obtain the highest safety coefficient, should not exceed 3/4”. 
It was also found that, contrary to the common view, cement concrete roads become 
slippery when wetted and their safety coefficient is not high. 


INTRODUCTION 


One of the main causes of vehicle accidents is skidding. In order to improve the traf- 
fic safety of the vehicle, it is thus important to know the degree of friction and 
its variation under given conditions. In 1959, measurements were carried out by 
the Division of Highway and Soil Engineering of the Technion with the aid of the 
Leroux “rugosimeter” (Figure 1), the principal part of which is a pendulum pro- 


Figure | 
Leroux rugosimeter. P—Pendulum. 1—Initial position. 2—Lowest position. 3—Rise. M—Centre 
of gravity of pendulum. a—Rubber lining. D—Length of contact. 
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vided with a shoe lined with standard-quality rubber. The pendulum is swung from 
a fixed height, losing part of its kinetic energy through friction at the road surface 
and as a result rising to a smaller height in the second phase of the swing. The 
percentage difference in height (H-h)/H, read directly from the scale, represents 
the proportion of energy absorbed. 

This method, unlike most others involving the determination of forces (measured 
either on the wheels of the vehicle, or on a special fifth wheel) is based on the deter- 
mination of energy absorbed by friction. 

In practice, the test consists in determining the percentages of absorbed energy 
for dry and wet surfaces. These are denoted by EF and e respectively, and referred 
to hereafter as the “‘wet’’ and “dry” roughness factors (coefficients de rugosité). 


The ratio K, = e/E, giving the decrease in friction due to wetting, is defined in 
the Leroux method as the ‘‘safety coefficient’”’ (coefficient de securité) of the road. 
The definition is based on the assumption that dry friction is in most cases high 
enough (according to Leroux, E = 85%) to eliminate the danger of skidding, even 
with worn tyres, and that the main source of this danger lies in the sharp decrease 
in friction following sudden wetting and in the inability of the average driver to 
estimate this decrease. The above assumption is based on Continental conditions, 
involving unexpected change in road surface characteristics due to rains during the 
warm season. 


Leroux’ classification of road safety according to the allowable driving speed is 
as follows: 


Allowable speed High Moderate Low Danger 
(km/hr) (100) (80-60) (60) 
K, 0.85 0.85-0.75 0.75-0.60 0.60 


The above figures correspond to the following values for e: 


Allowable speed High Moderate Low Danger 


e 70 70-60 65-50 50 


The relation between safe speeds for dry and wet roads (V and vy respectively) is 
given by 


4£mV2 =4 K, mv2 


which follows from the principle underlying the rugosimeter test. 
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In practice, the following average speeds were observed by Leroux on roads with 
Ki ==10.70: 


Dry road 120 110 100 90 
Speed (km/hr) 
Wet road 100 92 83 75 
EXPERIMENTAL 


In the present study, it was decided to base the classification of the road surfaces 
tested on the roughness factors proper as well as on the safety coefficient, for the 
following reasons: 


(a) A dry roughness factor of 85% was found only in 40% of the cases under 
consideration. It should be noted that measurements were mostly taken on main 
roads in a better state of repair than others. 


(b) A wet roughness factor of above 70% was found only in 16% of the cases. 
This means that even with a satisfactory safety coefficient, most of the roads under 
consideration were still unsafe for driving in wet weather. 


(c) Local climate differs from that of the Continent mainly in the existence of 
distinct wet and dry seasons, with no rains at all during the latter and little or no 
danger of unexpected changes in road surface characteristics. 


The following table shows the allowable speeds as determined on the basis of 
the above assumptions: 


E 85.0 85.0—80.0 80.0 

e 65.0 65.0—60.0 60.0 

Ks 0.75 0.75—0.70 0.70 
Allowable speed High Moderate Low 


The allowable speed category was determined as the lowest corresponding to the 
coefficients observed. 


A total of 26 points were tested (see Table J). In each case dry and wet measure- 
ments were taken in the direction of traffic and against it, and, at the same time, 
in the direction of the longitudinal slope of the road and against it. In the wetting 
tests a 5-10 mm layer of water was used. Results are summarised in Figure 2 giving 
the averages of the rugosimeter readings in its different positions. The test stretches 
were classified in groups according to the type of the road surface (Table I). 
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TABLE I 


Road 


Location 


Point Type of surface or trade name Age (months) General state 
1 Haifa—Nazareth Opposite “‘Zera’’ factory “Standard” 36 Cracked and patched 
Dp Haifa—Nazareth Near Acre Road intersection “Standard” 120 Widespread cracking 
3 Haifa—Acre 500m north of Acre Road intersection “Standard” 60 Cracked and patched 
4 Haifa—Hadera Near Geva Carmel village Bituminous concrete 12 Good 
5 Haifa—Hadera Near Atlit prison camp site “Coarse Shellmac”’ 12 Limited cracking 
6 Haifa—Hadera Near Faradis village, km. 67 Bituminous concrete 12 Good 
7 Haifa—Hadera Near Sitrine Agricultural School Bituminous concrete 12 Good 
8 Naveh Shaanan Upper stretch “Coarse Shellmac” 36 Cracked 
9 Haifa—Hadera 500m from Hadera Town “Standard” 36 Cracked 
Planning Boundary 
10 Hadera—Netania 500m from Hadera Road intersection Bituminous concrete 6 Good 
ft Hadera—Netania 200m from Giv’at Olga intersection Bituminous concrete 6 Good 
12 Hadera—Netania Near ‘“‘Abir’’ factory “Coarse Shellmac”’ (Over 3 years) Cracked 
13 Tel Aviv—Netania 100m from Wingate Institute “Coarse Shellmac” (Over 3 years) Cracked 
entrance 
14 Tel Aviv—Netania 200m from Tel Aviv Town Bituminous concrete 12 Good 
Planning Boundary 
15 Kiryat Bialik Residential street near “Ahava’”’ Cement concrete 240 Fair 
School 
16 Kiryat Bialik Residential street opp. ““Ahava” Cement concrete 240 Fair 
School 
17 Kiryat Bialik Residential street opp. ““Ahava” Sealing coat (over “Standard”) 36 Good 
School 
18 Kiryat Haim 300m from seashore Cement concrete 240 Fair 
19 Kiryat Haim 300m from seashore Sealing coat (over cement 36 Good, rough surface 
concrete) 
20 Beit Lid—Petah Tikva 500m from Tel Mond Road intersection “Coarse Shellmac”’ (Over 3 years) Fair, patched with bitumen 
21 Ramla—Jerusalem Near Pedayah village intersection “Standard” (Over 3 years) Irregular surface 
22 Ramla—Jerusalem 500m from Re’em village intersection “Coarse Shellmac” (Over 3 years) Widespread cracking 
23 Haifa—Nazareth Near Acre Road intersection “Standard” 4 New surface, still rough 
24 Haifa—Nazareth Near entrance to Sha’ar Ha’amakim Sealing coat 24 Limited cracking 
25 Haifa—Nazareth 200m from Tiv’on intersection “Fine Shellmac” 60 Cracked, patched with 
: bitumen 
26 Haifa—Nazareth “Fine Shellmac” 60 


Near entrance to Ramat Ishai village 


Fair, patched with bitumen 
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Figure 2 
Distribution of test points according to allowable speed 
TABLE IL 
Type of surface bey. Max. size of Bituminous Appr. percentage 
(or trade name) Description aggregate binder of binder 
1. “Standard” Dense graded (old) 3/4” ASG. 4% 
80/100 
2. Cement concrete iP? No bituminous coat 
3. Sealing coat over Dense graded Wie A.C. (1 kg/m?) 
cement concrete 80/100 
4. Bituminous concrete Dense graded il? A.C, 6% 
80/100 
5. ‘Coarse Shellmac”’ Open graded 2A 2a A.C, 4% 
80/100 
6. ‘Standard’ (2 weeks old) Dense graded 3/4” A.C, 4% 
80/100 
7. Sealing coat Dense graded lip A.C, (1 kg/m?) 
80/100 
8. “Fine Shellmac”’ Open graded ie A.C. 4% 
80/100 


DISCUSSION OF RESULTS (Table III, Figure 3) 

(1) The dry roughness factors found were fairly high, the average minimum 
observed being 77.5%, average maximum 93.4% and total average — 83.9%. The 
corresponding figures for e were 49.9, 74.6 and 63.3% and for K, — 0.635, 0.845 and 
0.75 respectively. 

(2) “Wet” readings showed higher dispersion than “‘dry” ones, which accounts 
for a correspondingly higher dispersion of the safety coefficient. 
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TABLE Ill 
Summary of results according to type of road surface 


a 


Point Grade 
Group __ Type of surface RS: E é K, = e/E (per cent) 
2 See. Be ae a ae a oe a ea ee 
1 “Standard”’ (old) 1 80.7 59.6 0.735 1.0 
2 VP 60.8 0.78 0.5 
9 82.9 58.4 0.70 — 
3 81.4 65.2 0.80 0.5 
21 85.0 65.7 0.775 4.0 
Average 81.4 61.9 0.758 
2 Cement concrete 15 81.4 57.6 0.71 — 
16 79.6 59.3 0.745 1.0 
18 85.7 54.4 0.635 0.5 
Average 82.2 Sy/al| 0.697 
3 Sealing coat over 17 82.4 70.1 0.845 “| RY 
cement concrete 19 82.1 68.0 0.825 0.5 
Average 82.3 69.1 0.835 
4 Bituminous concrete 4 80.4 69.2 0.87 1.0 
6 83.6 74.6 0.89 — 
a 86.5 62.3 0.72 
10 81.6 64.3 0.79 1.0 
11 83.3 66.9 0.80 ee 
14 82.5 63.9 0.77 — 
Average 83.0 66.9 0.81 
5 “Coarse Shellmac” 5 83.4 58.0 0.695 1.0 
8 78.7 65.2 0.825 4.0 
12 81.7 49.9 0.605 3.0 
20 87.3 58.9 0.665 5.0 
22 84.9 62.9 0.745 7.0 
13 85.2 54.5 0.640 0.5 
Average 83.5 58.5 : 0.696 
6 “Standard” (new) 23 89.3 72.0 0.805 — 
/ Sealing coat 24 90.8 57.6 0.64 5:5 
8 “Fine Shellmac”’ 25 93.0 65.6 0.705 7.0 
26 93.4 73.4 0.78 4.5 
Average 93.2 69.5 0.743 


(3) Comparison of the safety coefficients shows that high speed is permissible 
only on the new “Standard” surface (Group 6, i.e. 4°% of all cases), moderate speed — 
for groups 1, 3, 4, 8 (““Standard’’, sealed cement concrete, bituminous concrete and 
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100.0 


50.0 


Figure 3 


Distribution of types of road surfaces according to allowable speed 
Types of surfaces: 1. “‘Standard’’ (old). 2. Cement concrete. 3. Sealing coat over cement concrete. 
4. Bituminous concrete. 5. “Coarse Shellmac’’ (old). 6. “Standard” (new). 7. Sealing coat over 
Asphalt. 8. “Fine Shellmac’’. 


‘Fine Shellmac’’, a total of 35%), and only low speed for the remaining groups (over 
60%). This contradicts the common view that coating roads with a bituminous 
layer with coarse aggregate leads to increased friction and improves road safety 
during the rainy season. In fact, surfaces made with coarse aggregate (“Coarse Shell- 
mac’’) showed the lowest roughness factors and safety coefficients. It seems that the 
aggregate size providing the highest roughness factor and safety coefficient does 
not exceed 3/4’. 

(4) Cement concrete surfaces also showed low safety coefficients. This is also 
contrary to the common view that cement concrete roads are unaffected by wetting. 

(5) No difference was found between rugosimeter readings taken in the direction 
of the slope and against it. 

(6) Owing to the limited number of tests, no appreciable difference could be found 
between rugosimeter readings taken in the direction of the traffic and against it. 
This problem deserves further study. 


CONCLUSIONS 
The present study, the first road survey of its kind in Israel, has shown that local 
road surfaces are unsuitable for high-speed driving. It seems that the use of coarse 
aggregate does not contribute materially to driving safety, and it would be advisable, 
in future studies, to consider the influence of aggregate size and type of bitumen on 
the safety coefficient. 
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It is also advisable to determine (a) whether the decrease in roughness factor due 
to wetting is gradual, (b) the boundary between “‘wet” and “dry” in terms of driving 
safety, (c) the role of the texture of the road surface in this respect. 

Another problem deserving attention is the influence of the age of the surface layer. 
The tests reported here did not provide enough data on this problem, although it 
is already clear that it is not the chronological age alone that counts, but also the 
rate of wear of the road and its general condition. 


REFERENCES 


Leroux, M., 1957, La Route et la Securité—Etude de L’ Adhérence, Chartres. 

2. THOMANN, E. AND ZIpKEs, E., 1945, Die Reibungskennziffer und ihre Anwendung bei der 
Beurteilung von Strassenbelagen, Mitteilungen aus dem Inst. fuer Strassenbau der Eidg. Tech. 
Hochschule, Zurich. 

3. Grime, G. AND Gigs, C. G., 1954-5, The skid resisting properties of roads and tyres, Proc. 
Inst. Mech. Engrs., Automobile Division, (1) 19-30, Discussions 45-48. 

4. Moyer, R. H., 1934, Skidding characteristics of automobile tires on road surfaces and 
their relation to highway safety, Bulletin 120, Jowa Engineering Experiment Station, Iowa 
State College, 33, No. 10. 

5. 1938, Internationaler Staéndiger Verband der Strassen-Kongresse, VIII Kongress, Haag, 

Untersuchung und Messung der Beschaffenheit eines Strassenbelags unter dem Gesichtspunkt 

der Glatte und des Widerstands gegen Rutschen (No. 68-77). 


= 


MEASURES OF DEFORMATION IN THE STRAINED AND IN THE 
UNSTRAINED STATE* 


Z. KARNI AND M. REINER 


Technion—Israel Institute of Technology, Haifa 


ABSTRACT 


Expressions for the Green-measure of deformation in the strained state and of the Al- 
mansi-measure of deformation in the unstrained state are derived. 


1. Green and Almansi®’”’'**~ “postulated measures of finite deformation for which 


Murnaghan; derived tensor expressions. Let ,a be the ‘“‘initial’’ coordinates of 


a particle, x; its ‘“‘final’? coordinates in a Cartesian system, then the Green 
tensor 7°?:77~8 is 


Ox, Ox : (1.1) 
ny. ese Gi | (ee eed ae 
ina) 2 (= 0,a ij ) 


and the Almansi tensor is2 


€(x;) = 4 @ sue si) (1.2) 


Ox; Ox; 


As indicated by the placing of the subscript, the first refers to the unstrained state, 
the second to the strained state. This implies two difficulties. First, if it is desired 
to calculate the stress from a stress-strain relation, the final, strained state form must 
be used because only the final state is a state of equilibrium. Secondly, the initial, 
unstrained state of a body is the one given, and in certain problems of elasticity we 
must start from it. We therefore require expressions for y,;,; and for ,,e. These will 
presently be derived. 

2. Consider a spherical body of unit radius, which is subjected to a homogeneous 
deformation. Let the strain ellipsoid have the half axes 1,, /,, 1,, or in general I, 
then the principal strains are given by 


2e=1-1/P | (2.1) 
Or 


1—2¢ = 1//? (2.2) 
2n+1 =P 


* The research leading to this paper has been sponsored in part by the Air Research and Develop- 
ment Command, United States Air force under Contract AF 61(052)-223 through the European 


Office, ARDC. 
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We see that 
(1 — 2e)(2n +1) =1 (2.5) 


This equation can be generalized to a tensor expression by writing ¢;; for ¢, 1; for 
and 0;; for 1. We thus get 


(ig — 26;g) (2ngj + 9p;) = 9: (2.4) 


Now introduce from Eq. (1.2) 


gir Boggart a - (2.5) 
This gives 
(2ng; + 5p) we & = 5, (2.6) 
or 
pi a = =4(6, = aa 2) (2.7) 


On multiplying both sides of Eq. (2.7) by oh oe we obtain 
I n 


Narn aiass ta! aa 1mhd Cu moe ae 
x; OX, O,a O,a Oa 0,a Ox; Ox; Oa 0,a 


(2.8) | 


We now contract first on the indices i and k by putting i = k = A, and then onthe | 
indices ] and n by putting ] = n = o, bearing in mind that 


i Mir C5 Gksacn 
ax, a ox, a re 


and find (writing finally i for m and « for o) 


Ox; Ox; 
U u(5 a aa = su) (2.10) 


Similarly we find 


0;a O;a 
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3. We apply these considerations to the problem of simple shear given by 


NORE SP By SS I are (3.1) 
a=x-—TIy, b=y, c=z 
We find from Eq. (1.2) 
0 LB 0 | 
2; = 4 | i es (3.2) 
nee 0 0 | 
and from Eq. (1.1 
q. (1.1) is Romi 
ii = 4 if ie 0) (38) 
0 0 Oe 


The usual procedure for deriving n;; from the expression for ;;y is by carrying out 
a rotation through the angle a, which brings the initial into the final state+, where 
in accordance with Love! (1927) Art. 3 


tan geo | 2 (3.4) 
The rotation matrix is 
| cose sing 0 | 
Rj; | —sing cosa 0 | (3.2) 
0 0 0 
We find 
\ re oT OP 
i 
I OF =O =i 


This result can be obtained directly from the application of Eq. (2.10). While Love’s 
method is applicable in this simple case, Eq. (2.10) only can be used in other more 
complicated—for instance—three dimensional problems. 

Finally we find from Eq. (2.11) 


| =|—— IL 0 
ye = 4 | tO O | (3.7) 
tO OL 0 | 


4, Murnaghan and, following him, Michal have named the Almansi tensor *‘Eu- 
lerian’’, and the Green tensor ‘‘Lagrangean’’. They have chosen this terminology 
in analogy with the practice in hydrodynamics. However, while in hydrodynamics 
these two terms result from two different viewpoints of the same thing, the Green 
and the Almansi tensor define two different things as is obvious when we compare 
Eqs. (3.2) with (3.6), or (3.3) with (3.7). This becomes strikingly evident when we 
calculate from the classical stress-strain relations the stresses connected with the 
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simple shearing deformation. In this case the use of the Green measure requires 
the existence of a tension in the direction of the shearing displacement (x), while 
the use of the Almansi measure is legitimate only if there is a pressure in the direction 
of the shearing gradient (y). This is not always realized. For instance, Sneddon 
and Berry> (1958) speak of ‘‘the components of the strain tensor in terms of the 
components of the displacement vector (uv) with respect to the (original state) S’’, 
and ‘‘in terms of the components of wu with respect to (final state) S”’, as if it were 
the same strain tensor. 
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ON THE DETERMINATION OF AIRY POLYNOMIAL STRESS FUNCTIONS 


Zv1 HASHIN 
Division of Mechanics, Technion-Israel Institute of Technology, Haifa* 


ABSTRACT 


A straightforward method is presented for the determination of polynomial Airy stress 
functions, which give solutions for continuously loaded rectangular strips. This is an ex- 
tension of a method given by Ching-Yuan Neou in a previous paper. 


1, INTRODUCTION 
The two dimensional problem of the second kind of the theory of elasticity has been 
reduced by Airy to the following boundary value problem: Solve the equation 
04F O4F O4F 
Are ee ey 1.1 
VF aNes 2 ayn Oe 0 (1.1) 
so that on the boundary of the region the stress vectors assume specified values. 
The stresses are defined by the second derivatives of the stress function 


02F 
6, = by? (1:2) 
02F 
%, = (1.3) 
02F 
= 1.4 
“y dxdy ae 


In the standard texts on the theory of elasticity’’* we find so-called ‘‘elementary”’ 
solutions for the two dimensional problem. These are Airy stress functions which 
are expressed by polynomials in x and y, which in general give solutions for rectan- 
gular domains. 

The chief value of these solutions lies in their simplicity, but the usual methods 
for obtaining them are not straightforward. The solution is either guessed or found 
by combining known solutions so as to satisfy the boundary conditions. 

In a recent paper by Ching-Yuan Neou3, a direct method for determining such 
stress functions is described. The method still remains somewhat arbitrary in that 
no specific reason is given for the choice of the degree of the polynomial stress 


* Now at Division of Engineering and Applied Physics, Harvard University, Cambridge, Mas- 
sachusetts, U.S.A. 
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function. This has also been pointed out in the discussion on the above mentioned 
paper‘. 

In the following this point will be cleared for the case of rectangular domains. 
We also believe that the method here proposed, for fitting the polynomial to the 
boundary conditions, reduces substantially the amount of computations required 
by Neou’s method. 


2. METHOD OF SOLUTION 


We shall here deal with problems of the second kind; in these the loading on the 
boundary is prescribed. Instead of imposing boundary conditions on the stresses, 
which amounts to prescribing the second derivatives of F on the boundary, we shall 
use a different kind of boundary conditions which will be found very convenient 
in our case. 


Figure 1 


Let AB be part of a boundary curve of a two dimensional region (Figure 1). The 
length of the curve s is measured from A to B, counter clockwise. Let x,, y, be the 
coordinates of any point B which is at distance s from the reference point A. Let 
x,y be the coordinates of an arbitrary point on the boundary between A and B 
and T,, T, the components of the boundary stress vector at this point. It may 
be proved'then that the stress resultants S, and S, and the moment M are given by 


s 


A. OF 
Si7 = | T,ds = —— 
ea (2.1) 
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The values of F and its derivatives are thus determined at the arbitrary point B 
assuming that their values at A are zero. This last assumption can always be made as 
it is equivalent to the superposition of an arbitrary linear function on F. This will 
not change the stresses as they are expressed by second derivatives of F. 


The signs of S,, S, and M are determined as follows: the positive sense of traversing 
AB is counter-clockwise; a force is positive when its arrow points outwards, a moment 
is positive when it turns counter-clockwise. S, and S, will be positive when they 
point outwards while the contour is traversed in the positive sense, or when they 
point inwards while the contour is traversed in the negative sense. The same rea- 
soning applies to the moment; it will be positive when the sense of traversing AB 
and its sense of rotation coincide. 


We write the stress function F(x,y) in the form 


@ 


Fey) = ET Caety" (2.4) 


m=o h=o 


Introducing this into (1.1) we obtain, after equating to zero coefficients of equal 
powers, the following systems of equations by which the C,,, are interrelated.* 


(m +2) (m+1) m (m—1) Cysa2n-2 + 2m(m—1) n (n—1) C,, + 
(n+2) (n+1)n (n—1) C,,-2.n+2 = 0- (2.5) 


We now arrange the coefficients C,,,, in an infinite rectangular matrix 
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It is easily perceived that any three coefficients which are related by an equation 
of type (2.5) are found in the matrix by going along a diagonal from left to right 
upwards. For example, C4,, C24, Cog which are encircled in the matrix. 


Let us assume that the polynomial (2.4) does not contain powers in x higher 
than k. Then all coefficients in the rows for which m > k will vanish. For k = 5, 
say, all coefficients below the horizontal line indicated in the matrix are set equal 
to zero. Applying equations (2.5), a little reflection will show that only the coefficients 
inside the region indicated by the double line can be different from zero. 

In general, when the x powers are terminated at k, this region may be constructed 
as follows: in the kth line go 4 spaces horizontally to the right, then 2 spaces up- 
wards, then two spaces to the right etc., till the row m = 0 is reached. We thus ob- 
tain a step region as shown above. When k is even it will be found that the height 
of the last step is only one. 

In order to determine completely the stress function we shall have to determine 
all coefficients within the region and we shall accordingly need as many equations. 
It is easily proved that when m is terminated at k the number of the C,,,, to be de- 
termined is 


R 
R 


M(k2 + 8k +7) for k odd (2.6) 
4(k2 + 8k + 8) for k even (2.7) 


The equations from which these coefficients are determined are supplied from 
two sources 

a. Compatibility relations (2.5) 

b. Boundary conditions. 

We shall first determine how many equations of type (2.5) may be obtained from 


a typical step region when m is determined at k. Again we have to distinguish be- 
tween the case k-odd and k-even. Two such step regions are drawn below 
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By proceeding from one diagonal to another it is found that the number of com- 
patibility relations which may be set up, by using the coefficients of a region, is 
essentially the sum of an arithmetic series. 

It is not difficult to prove that this number is given by 

Q = 4(k2 — 1) for k odd (2.8) 
QO =4k2 for k even (2.9) 

The number of equations which are still needed and have to be obtained from the 

boundary conditions is thus 
B=R-@Q 
In both cases it is found that 


B= 4k +1) (2.10) 


We now turn to the case of a two dimensional elastic body of rectangular form 
(Figure 2). 


Figure 2 


Let us assume that the faces y = +a are loaded by a continuous stress distribution 
giving a moment with respect to some fixed point on the boundary which is a poly- 
nomial of degree k in x. This is equivalent to assuming a stress o, in form of a poly- 
nomial of degree k — 2. Furthermore, it is assumed that the horizontal force S, 
is a polynomial of degree not higher than k in x. This by (2.1) (2.2) (2.3) is equiva- 
lent to prescribing F and 0F/dy on the boundaries y = +a. (0F/0x is already deter- 
mined by prescribing F on the lines y = const.) 

On the vertical sides x = +I the stresses shall not be specified. It will be demanded 
that these have resultants which satisfy the equilibrium conditions of the body. 
Thus our solution will only hold within the limits of validity of Saint Venant’s 


principle. 
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It is important to realize that specification of the resultants of the stresses on the 
vertical sides is not equivalent to additional boundary conditions, on F. When 
advancing around the contour along AB, through BC, say to CD—the passage of 
BC will add the moment of the stresses on it directly to the moment along CD so 
that we are left with conditions on F on AB and CD only. The same applies to the 
force resultants, i.e. derivatives of F. 


This very much simplifies the calculations. In usual methods of solution, when the 
boundary conditions are expressed in terms of stresses, these involve constants 
which are determined from integration over the vertical sides. This in general leads 
to cumbersome computations. 

F(x,y) as a polynomial of degree k in x, may be written in the form 


m=k n=l 


F(x,y) = Tis Y > Ce y" (2.11) 


m=0 n=0 


where the highest power of y is determined from the step region. Let the boundary 
value of F on y = + a be given by 


F(x, a) = > M,x' (2.12) 
on y = —a by: is 

F(x, —a) = So N,x" (2.13) 
Let S, = OF /dy on the boundary be given by - 

= ae i Lx" (y = +a) (2.14) 

0 r=k 

nD Pe (y = -a) (2.15) 
Differentiating (2.11) we obtain 

m=k n=l 
a7 2a, NC omX"Y" ? (2.16) | 


Equating (2.11) for y = +a to (2.12) we obtain k + 1 equations 


Y Cra” = My raze Oy, 2h. i (2.17) | 


In an exactly analogous way we obtain k + 1 equations for each of the remaining 


three boundary conditions, in all.4 (k + 1) equations. This is exactly the number of | 


equations still needed in order to determine the coefficients C,,,, of F(x,y). 
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We have thus arrived at the following conclusion: if on a rectangular domain 
the stresses are given and expressed as polynomials on two opposite sides, while 
the stress resultants are prescribed on the two other sides—the problem may be 
solved by an Airy stress function, in the form of a polynomial and the method of 
solution is that described above. 


It is evident that exact boundary conditions on stresses apart from those on two 
opposite sides, cannot be satisfied by our method of solution. These would make 
the number of equations greater than the number of unknown coefficients. 


3. SOLVED EXAMPLE: BENDING OF BEAM UNDER UNIFORM LOAD 


We shall now illustrate the method previously described by working out in detail 
the solution for a beam on two supports, uniformly loaded on its upper face. We 
shall assume that the beam is of unit thickness. The boundary conditions for F that 
are to be fulfilled are found as follows: We choose point A as zero point and traverse 
the contour clockwise. Because of the symmetry with respect to the y axis it will 
suffice to consider one-half of the beam. 


Figure 3 


On the vertical boundaries BC and EF the sum of the shearing stresses must 
equal the reactions p/ and the bending moments are zero. The boundary conditions 
on the horizontal faces are: 


On AB: y=a 


MG) =k). ¥ ea (3.1) 


LS ae, (3.2) 


S,(x) = Oy. 
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on CD: y=-a 


p? p? 
M(x) = F(x, -a) = —=- + M, = ——= (3.3) 
OF r 


As the moment in (3.1) is of the second order in x, we choose a stress function of 
order 2 in x. The step region of coefficients is then as shown below: 


Because of the symmetry with respect to the y axis, coefficients of odd powers 
of x may be equated to zero, i.e. 


Cho Cit = Gin = Ci, = 0 
Our stress function thus becomes 


F(x,y) = Coo + Cory + Cory? + Cosy? + Co4y4 + Cosy® (3.5) 
+ X2(Co9 + Cayy + Cooy? + Co3y3) 


OF 
aay. = Cor + 2Cory + 3Co3y2 + 4Cogy3 + SCosy* (3.6) 


+ x2(C21 + 2C22y + 3C23y?) 


From boundary condition (3.1) we obtain the following equations 


Coo + Cora + Cora? + Co3a3 + Coyat + Cosa? = 0 (3.7) 
Coo + Coa + Cy2a? + C303 =-§ (3.8) 
from (3.2) 

Coo — Co14 ar Co2a2 —s C343 ap Co4a4 = Cos5a° = (3.9) 
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Crp — C213a + Cz2a2 — C,,3a3 = 0 (3.10) 
from (3.3) j 
Co, + 2Co2a + 3Co3a2 + 4Co4a3 + 5Cy5a4 = 0 (3.11) 
Cy, + 2C,2,a + 3C,3a2 = 0 (3.12) 
from (3.4) 

Co, — 2Co2a + 3Co3a2 — 4Co4a3 + S5SCo5a4 = =0 (3.13) 
C2, — 2C,2,a + 3C,3a2 Sper ib) (3.14) 


We have thus obtained 8 equations for the unknown coefficients. Two addi- 
tional equations are obtained from the compatibility relations (2.5). This checks 
with (2.9) which yields Q = 2 for k = 2. These 2 equations are: 


Cr, + Soa => 0 (3.15) 
Co3 +°5Co5:= 0 (3.16) 
Equations (3.7) — (3.16) are easy to solve, especially when suitable pairs are 
added and subtracted. The unknown coefficients are: 
pl? fs P 
Coo F Tiss Cos Sadaeaar 
3140011 i -wep 
Co1 =P (aa) Cr Aid 
3p 
Cor. = 0 . Cy. = - 8a 
Fn P y fe 
Co3 a 8a3 C22 = 9 
P 
Cos = 0 D3 ae 8a3 


Of these the first two are not of interest for the determination of stresses. The 
stress function is accordingly: 


@ ealeulaPrny WIG) Ho aud, iD, 8 Pe Pag 7 
Fe = (ah Ey eae er re mei ey 


The stresses are now easily found by differentiation. We thus obtain 


2) By oe 3.18) 
ail Vesa 4 BEG @2 ( 

_ 1 ae ee ae 3.19 
= i (= Ne 2) (3.19) 


Clo p.41) 
eee. er (: i a) (3.20) 


a2 
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4. CONCLUSION 


We have shown how a certain class of boundary value problems of the two di- 
mensional theory of elasticity, can be solved by means of polynomial stress functions. 
The solution can be constructed easily and in a straightforward way in all cases 
where the domain is rectangular and continuous stress distributions—in terms of 
polynomials—are prescribed on two opposite sides and stress resultants on the two 
others. We thus conclude that the method of solution is of practical value only 
when the ratio between length and breadth of the domain is sufficiently large to 
permit the application of Saint Venant’s principle. 

When the stress distribution on the long sides is piecewise continuous, a solution 
can be constructed by cutting the domain vertically at the point of discontinuity 
and solving for each part separately. It will then be demanded that the stress re- 
sultants on the cut be equal to the statical values determined from equilibrium of 
the parts. It must be added that the practical value of such a solution is questionable 
as it will take no account of the stress concentrations that will usually appear at the 
point of discontinuity. For this kind of boundary value problem a better solution 
would be obtained by the well known method of Fourier Series. 
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RESISTIVITY CALCULATIONS FOR SEMICONDUCTOR CRYSTALS 
GROWN BY THE CZOCHRALSKI TECHNIQUE* 


B. PRATT 
Department of Metallurgy, Technion—Israel Institute of Technology, Haifa, 


M. Y. BEN-SIRA 
Department of Physics, Technion—Israel Institute of Technology, Haifa 


ABSTRACT 


The derivation of a practical formula for the calculation of the top resistivity of semi- 
conductor single crystals grown by the Czochralski technique is presented. This formula 
contains measurable quantities only, e.g. mass, resistivity, mobility and segregation co- 
efficients. It is shown that, in practice, aiming at the top resistivity is enough in order to 
bring the crystal into a required resistivity range. 


I. INTRODUCTION 


Silicon and germanium are most commonly used in the manufacture of semicon- 
ductor devices. Single crystals of either n-type or p-type material are used for these 
devices. These crystals are of a wide range of resistivity differing from each other 
according to the resistivity requirements of the various devices. A widely used 
technique for producing single crystals with the required type and range of resistivity 
is the growing of crystals by the Czochralski technique!. Here the crystal is pulled 
out of a melt of the starting material doped by the right amount of appropriate im- 
purities. The calculation of the amount of “dope” needed involves rather compli- 
cated computations. 


In this paper a simple and practical formula is presented, useful for the dope cal- 
culation of crystals grown by the Czochralski technique. The formula contains 
measurable quantities only, e.g. mass of melt, resistivity, mobility and segregation 
coefficients. 


* This work was carried out during the stay of the authors with the International Rectifier Cor- 
poration, El Segundo, California. 
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2. THE CZOCHRALSKI TECHNIQUE 


In the Czochralski method, the crystal is pulled out of the molten material with 
the aid of a “‘seed’’. The seed is a piece of single crystal of the same material as that 
of the melt and is held to a special shaft so that its desired orientation will be parallel 
to the surface of the melt. The melt, contained in a special crucible, is kept slightly 
above the melting point. The seed is dipped into the melt and slowly pulled out, 
with the solidified material, at a speed determined by the rate of solidification and 
the desired diameter of the grown crystal. The melt is kept in an inert gas atmosphere 
or under high vacuum, in order to avoid oxidation at the high melting 
temperature. The melt and the seed are both rotated while pulling the crystal, to 
maintain uniform heat distribution and crystal growth. These conditions make it 
possible to grow a single crystal as the solidification takes place in a desired orient- 
ation. 


The melt contains a certain concentration of impurities, usually added according 
to the starting material and the required crystal resistivity. When the material so- 
lidifies the concentration of impurities in the solid will be only a certain quotient K 
(the segregation coefficient) of the concentration of impurities in the melt, depend- 
ing on the semiconductor and the impurity. Thus while pulling the crystal, the im- 
purity concentration in the melt will change continuously and with it the impurity 
concentration in the crystal will also be changing. Since the resistivity of a semi- 
conductor crystal is inversely proportional to the impurity concentration, the 
resistivity will vary along the crystal, forming a “resistivity profile’. This profile 
can be computed with th aid of the concentration equation?: 


Ons ke —g)*! 


where c is the concentration of impurity at a point of the crystal, k is the segregation 
coefficient of the impurity in the semiconductor, c, is the concentration of impurity 
in the melt before crystal growth starts and g is the fraction of melt solidified up 
to measured point on crystal. This rather complicated formula requires, in practice, 
the use of special tables and calculations in order to compute the amount of im- 
purities required for a certain profile. 


From this formula it is also seen that fork < 1 the concentration c will 
increase with g; and for k > 1, c will decrease with increasing g. Actually, for semi- 
conductors such as silicon and germanium and for the commonly used doping 
materials, phosphorus, arsenic and antimony (n-type impurities) and boron, alu- 
minium and indium (p-type impurities), the segregation coefficient is smaller than 
unity. (The case of boron impurity in germanium is an exception, the segregation 
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coefficient being ~20). This means that when there is a large excess of one type of 


impurity only, the top of the crystal will be the highest in resistivity and the bottom 
will be lowest. 


In practice, the impurities are added into the melt in a piece of semiconductor 
with known very low resistivity, cut from a specially grown crystal heavily doped 
with a required kind of impurity. This way the amount of impurities added is easily 
control ed. For carrying out the computations of the resistivity profile with the 
aid of the above mentioned formula, the number of impurities either per cm3 or 
per 109 atoms of the semiconductor should be known. This quantity can be derived 
from the measurable resistivity of the semiconductor and the dope. 


3. CRYSTAL PROFILES AND RESISTIVITY RANGES 


One of the determining factors in the making of semiconductor devices is the re- 
sistivity of the starting material, for which the requirements vary according to the 
device. However, quite a wide range of resistivity is allowed for each particular 
case. The change of resistivity between the limits of these resistivity ranges is 
similar to the change of resistivity along a crystal grown by the Czochralski 
technique. In general, most of a crystal with top resistivity equal to one of the limits 
of a required resistivity range, will fall into that range. Therefore, if the doping of 
the melt is carried out so as to bring the top resistivity of the crystal to suit the 
value of the appropriate limit of the resistivity range, the rest of the crystal will 
then fit into the resistivity range in question. Thus for most practical cases it is 
enough to consider only the top resistivity of a crystal. 


A useful and practical formula is hereby derived for the calculation of the top 
resistivity of grown crystals. The formula contains measurable quantities of the 
used materials. 


4. DERIVATION OF THE FORMULA 


The conductivity o of a semiconductor is given by: 
(1) o = ent = -; orn = — 


where e is the charge of the electron, n is the number of carriers per unit volume 
(thus practically the number of excess impurities per unit volume), 4 is their mobility 
and r = 7 is the resistivity of the semiconductor. A semiconductor of weight W 
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and density d, which contains only one impurity element — associated with shallow 
levels, the sample in question being in its extrinsic range — will then have N car- 
riers given by: 


(2) N =—— (W/d being its volume). 


If a melt contains several, say m, such components the total number of carri- 
ers is: 


. Nee Sony ay ut 


> 
i=1 i=1 Cary; 


and here impurities of opposite type should have opposite sign in the summation, 
carried along with the electric charge. 


Only a fraction of the carriers from each component will enter from the melt 
into the solid crystal, namely K; x N; where K; are the respective segregation co- 
efficients. If it is assumed that the respective initial concentrations of impurities are 
maintained throughout the pulling process, then the resistivity along the crystal 
would stay constant from top to bottom and the total number of carriers in the 
crystal would be, by (2) and (3): 


(4) gia Ne cs Sokol 


edrsis ,<1 | edr,pt one hedge ipt S stands for the solid crystal). 


However, during the process of pulling the respective concentrations change con- 
tinuously and the resistivity r, in formula (4) will represent only that of the top of 
the crystal. This is the point at which the conditions postulated in the derivation of 
formula (4) are valid. 


Assuming the same density, at the melting point, for both solid and molten 
silicon (the actual differences being very small) and reducing by the factor ed on 
both sides of (4) the following practical formula is derived: 


(5) sergiilivami Ms 


i=1 nace rsUs 
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Here the summation is carried out over all the components of the melt, including 
the doping material. The n- and p-type components should be given opposite signs 
in the summation. As a rule it is more convenient to give the plus sign to the type 
of the required crystal. The fact that mobility is changing with the semiconductor 
resistivity and that it is different for n- and p-type components should be taken 
into account when using formula (5). These values of the mobilities are readily found 
in literature; and handbooks4. 


5. APPLICATIONS 


In practice, satisfactory results were achieved by applying formula (5) to doping 
calculations for silicon crystal growing, where very low-resistivity silicon pieces 
were used as doping material. A common case is that of doping a silicon melt of 
known resistivity and kind of impurity with the aid of such very low-resistivity 
silicon of known resistivity and kind of impurity. In this case, formula (5) becomes 


WN, We 
lulu Tplp rss 


(6) Ky 


where subscripts M, D and S stand for melt, dope and solid respectively. The only 
unknown in the formula in this case will be, most likely, the weight of the doping 
material. The same holds in cases where the melt is composed of more than one 
ingredient and then the resistivity and kind of impurity of each weighed component 
should be carefully evaluated. Here again the amount of known doping material can 
be calculated for a required top resistivity of grown crystal. 


Formula (5) can also be applied in the case of blending of low resistivity p-type 
and n-type silicon of known resistivity and kind of impurity in order to get a certain 
top resistivity of the crystal. Here either the weight of the p-type silicon or that of 
the n-type silicon will be the unknown. This blending is permissible only where 
highly compensated silicon is allowed. 


It should be noted that in the case of a very highly compensated component in 
the melt, results of resistivity calculations by formula (5) will not be accurate because 
the respective W/ru term will represent only the excess impurity situation and not 
the actual impurity contribution of that particular component to the melt. 


ACKNOWLEDGMENT 


The authors wish to thank Mr. L. S. Pelfrey of the International Rectifier Corpo- 
ration for useful ideas and helpful discussions on various aspects of this work. 


108 


> wowN 


B. PRATT AND M. Y. BEN-SIRA Bull, Res. Counc, of Israel 


REFERENCES 
Hunter, L. P., 1956, Handbook of Semiconductor Electronics, McGraw-Hill, pp. 6-15—6-31. 
PFANN, W. G., 1958, Zone Melting, Wiley, p. 10. 
DesyE, P. P. AND Konane, T., 1954, Phys. Rey., 94, 724. 
Hunter, L. P., 1956, Handbook of Semiconductor Electronics, McGraw-Hill, p. 7-5. 


QUANTITATIVE DETERMINATION OF FREE SULPHUR IN HYDRO- 
CARBONS BY THE BENZOIN METHOD 


A. B. STERN AND M. BABITz 


Chemical Testing Laboratory, Technion-Israe Institute of Technology, Haifa 


ABSTRACT 


A quantitative test for free sulphur in hydrocarbons is based on conversion of the free 
sulphur into H,S, determined colorimetrically as methylene blue. The reaction is sen- 
sitive to 0.5 y sulphur and unaffected by sulphur compounds normally occurring in pet- 
roleum. The method is especially suited for free sulphur determination in coloured pet- 
roleum products. 


INTRODUCTION 


Quantitative determination of free sulphur in oil products is of considerable im- 
portance in the treatment processes necessitated by its deleterious properties (such 
as corrosivity, lead susceptibility, etc.). Various authors have attempted to develop 
suitable methods, some of them by adapting qualitative tests involving the use of 
metals. For example, Garner and Evans! used copper bronze, and Stehr and Levin? — 
copper gauze, while Claxton and French? attempted to adapt the A.S.T.M. copper 
test+ for quantitative purposes. Others* resorted to mercury, which is a highly 
sensitive qualitative reagent. 


In all these methods (based on the conversion of metals into sulphides, reduction 
to H,S and determination of the latter), the analytical procedure is extremely com- 
plicated and time-consuming. It was found that peroxides have an inhibitory effect 
on copper and react with mercury; they also interfere in the reversed plumbite re- 
action suggested by Wirth and Strong’. In addition, these methods give poor results 
for low concentrations. 


The colorimetric method of Sommer® and Mapstone’s modification’ are based 
on the addition of pyridine and NaOH with sulphur precipitated in colloidal form, 
and invoive considerable errors in the micro range. A more accurate colorimetric 
method developed by Bartlett and Skoog!® is based on a reaction between free 
sulphur and NaCN in acetone, with NaCNS determined colorimetrically or spectro- 
photometrically (at 465 my wavelength) with the aid of FeCl,;. The method was 
found extremely accurate, provided dark containers are used and readings taken 
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instantaneously. However, in the case of red-coloured gasoline this method is liable 
to yield misleading results. 


Recently, polarographic methods were developed for the determination of free 
sulphur and sulphur compounds. These methods are most accurate, but necessitate 
expensive equipment. A rapid qualitative micro method suggested by Feigl'!is based 
on the conversion of elementary sulphur into H,S by oxidation of benzoin: 


C;H<CH(OH)COC,H; + 8 > CsH;COCOC,H. + HS 


The reaction according to Feigl takes place in the temperature range between 
140—150°C (and is extremely rapid, with H,S detected by means of lead-acetate 
paper. The authors attempted to adapt the colorimetric method suggested by Budd 
and Bewick!? and Jacobs et al.13 for the determination of H,S obtained by the 
Feigl method. 


The objects of the present study were: 

(a) ascertaining whether Feigl’s reaction is quantitative; 

(b) determination of optimum conditions; 

(c) development of an absorption technique permitting quantitative determina- 
tion of the H,S; 


(d) designing a simple apparatus for this purpose. 


REAGENTS AND APPARATUS 


1) Cd(OH), solution for absorption of H,S. 


2) Amino acid and FeCl; solutions for detecting H,S by the methylene blue 
method!3. 


3) (NH4)2HPO, solution, 400g in 805 ml water. 

4) Concentrated sulphur solution, 1 g sulphur per litre benzene (‘‘Analar’’). 
5) Standard sulphur solutions (200, 80, 40, 20, 10 and 5 y/ml) in n-heptane. 
6) Pure nitrogen, 99.9% (compressed in cylinder). 

7) Spectrophotometer—Beckman, Model B. 

8) Colorimeter—Lovibond tintometer (IP 17/52). 


9) Testing apparatus, as shown in Figure 1. 
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NITROGEN 


REACTOR 


GLYCEROL | 


BATH [" 


ABSORBER 


Figure 1 


PROCEDURE 


4g benzoin and 1 ml sulphur solution are placedin the reactor tube and 40 ml Cd(OH), 
solution in the absorption tube. (If room temperature exceeds 20°C the solution 
should be cooled). The tubes are joined and the reactor is connected with the nitrogen 
cylinder and immersed in a glycerol bath (initial temperature not to exceed 100°C). 
Nitrogen is passed through at a rate of 5 bubbles per second and the bath heated 
gradually to 180°C, after which the nitrogen stream is maintained for 5 more minutes. 
The apparatus is then dismantled and 0.6 ml amino acid solution are introduced 
into the absorption tube through the gas inlet, followed immediately by a drop of 
FeCl, solution; the tube is stirred and allowed to stand for 2 minutes, after which 
1 ml ammonium phosphate is added, the contents transferred to a 50 ml volu- 
metric flask and filled to the mark with distilled water. After 5 minutes the colour 
is measured in a spectrophotometer (670 my, 1 cm cell) or a Lovibond tintometer 
(1”, 4” and 4” cells, acc. to concentration), calibrated with standard solutions of 


Na,S. 


EXPERIMENTAL 


First tests were carried out at 150°C, as suggested by Feigl. Since at that temperature 
the H,S evolved is in danger of oxidation, a nitrogen atmosphere was provided 
between the reaction and absorption vessels. A pure (99.9%) nitrogen atmosphere 
gave satisfactory results, but when using commercial (97%) nitrogen, results were 
unsatisfactory even when all traces of oxygen had been removed with an alkaline 
solution of pyrogallol. 

The following conditions were examined: (a) temperature of conversion, (b) vol- 
ume and concentration of sulphur solution, (c) influence of various organic sulphur 
compounds. 
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Optimum temperature of the benzoin test 


The determination of free sulphur by the benzoin method was carried out at dif- 
ferent temperatures from 140°C to 190°C. Heating time varied from 30 to 40 minutes, 
from the moment of placing the reactor in the bath until the maximum temperaure 
was reached. It was found that optimum temperature was 180°C and that at lower 
temperatures the conversion was incomplete, while above 180°C no improvement 
was observed (See Table I). 


TABLE I 


Effect of temperature on sulphur conversion (without pre-heating of nitrogen) 


Final tem- 
Standard solution perature Lovibond reading Spectrophotometer 
(bath) reading 

S Volume (°C) S found Conversion S found Conversion 
(y/ml) (ml) (Y) (%) (y) (%) 

10 INcc: 140 5 50 4 40 

2 2 150 6 60 4 40 

He Ey 160 6 60 5 50 

% ? 170 9 90 6 60 

= is 180 10 100 7 70 

3 % 190 10 100 7 70 


Accuracy range 


Tests were carried out for determining the free sulphur in amounts of 2 to 80y in 
standard benzene solution, the volume varying from 0.2 to 4 ml. The heating time 
was again about 30 mins. and the temperature about 180°C. It was found that with 
concentrations exceeding 2y/ml the sulphur is capable of determination with an 
accuracy of 80-90% by tintometer and 70 to 80% by spectrophotometer. For 
lower concentrations a larger volume should be used. The maximum volume for 
the apparatus in question is 5 ml and the minimum detectable concentration 0.5 


y/ml. For concentrations exceeding 100 y/ml the solution should be diluted (see 
Table IT). 
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TABLE If 


Benzoin test on benzene sulphur solutions at various concentrations 


Final tem- 
Standard solution perature Lovibond reading ee 
(bath) 
S Volume (°C) S found Conversion SS found Conversion 
(y/ml) (ml) (y) (%) (Y) (%) 

2 0.2 180 2 100 1 50 

5 0.5 180 5 100 4 80 
10 1 190 9 90 7 70 
12 0.6 180 10 83 8.5 70 
20 1.0 180 18 90 15 75 
20 2 185 18 90 13 65 
30 3 180 27 90 24 80 
40 4 180 40 100 36 90 
80 1 182 70 87.5 60 75 
20 2 175 19 95 — eae 
20 2 180 19.5 97.5 = = 


Influence of various organic sulphur compounds 


Various sulphur compounds, such as mercaptans, sulphides, disulphides and carbon 
disulphide were tested. 


These compounds showed none or very little influence on the test. Olefinic sulphur 
compounds, apparently due to their chemical structure, are more capable of reaction 
and decomposition, and may give somewhat higher results due to liberation of free 
sulphur. Experiments also showed that free sulphur added to the above compounds 
could be detected quantitatively (Tables HI and IV). 


we 


Benzoin test on gasolines 


Different gasolines were tested both directly and after mercury treatment, in order 
to establish whether the sulphur detected is not due to decomposition of organic 
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TABLE III 
Benzoin test on mercaptans with and without free sulphur (Lovibond) 


(Temp, 180°C) 


Standard solution 


Free S N-buty! Volume S found 


(v) mercaptan (ml) (v) 
calc. as S (7) 


_- 12 1 0 

— 24 2 0) 

— 120 zZ 0 

— 600 I 0 

= 1200 2 ( 

10 12 2, 9 | 
20 24 2 18 
20 120 2 18 


TABLE IV 
Benzoin test on sulphides and disulphides with and without free sulphur (Lovibond) 


Compound Compound Free S Volume of 


calculated added solution S found Temperature 
as S (7) (v) (ml) (Y) (CC) 
Diphenyl | 
sulphide 160 = 1 0 175 | 
Diphenyl | 
sulphide 160 20 3 19 180 
Dibutyl 
disulphide 24 —_ p2 0 150 
Dibutyl 
disulphide 1200 — 2 <<) 180 
Dibutyl 
disulphide 300 20 25 19 180 
Diallyl 
disulphide 60 — 1 4 150 
Diallyl 
disulphide 60 — 1 10 200 
Diallyl F 
disulphide 600 — 1 24 180* 
CS» — — 1 0 180 


* Mercury treatment 
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sulphur compounds during the reaction. In straight-run and synthetic gasolines no 
free sulphur could be detected after treatment while in cracked gasolines it reap- 
peared in the amount of 30% of the previous result. 

These results are in agreement with the findings of the preceding section. 

Tests were also carried out on gasolines with admixtures of free sulphur. Results 
are given in Table V. 


TABLE V 
Benzoin test on various gasolines (Lovibond) 


S found 
Gasoli Volume S found S found after Hg 
eg (ml) (”) (mg/1) treatment Remarks 
(mg/1) 
79 octane 2 40 20 ot Cracked 
83 octane 
cracked 2 29 14.5 8.5 
83 2 27 IES) — Marketed motor gasoline 
83 ® 28 14 = 
83 octane 
+ 10yS D 40 20 9.5 
70 octane 2 0 0 — Straight-run 
70 octane 
+ 20S 2 18 9 0 
83 octane | 20 20 3355) Cracked 
100 octane l 0 0 0 Aviation gasoline, Sample A 
a 5 3 0.6 0 43 a eae A. 
| 2) 2 0 ia is aan 
ASTM deg. of corrosivity 2e 
BY 20 4 0 Aviation gasoline, sample B 
ASTM deg. of corrosivity 2e 
He 2 8 4 0 Aviation gasoline, sample B 
ASTM deg. of corrosivity 2e 
He 22 0 0 0 Aviation gasoline, sample C 
CONCLUSIONS 


Free sulphur in hydrocarbon solution can be determined quantitatively by benzoin 
conversion into H,S and colorimetric measurement of the latter as methyleneblue. 
The optimum conversion temperature is 180°C (measured in the bath). 

The minimum amount of sulphur detectable by this method is 0.5y and the volume 
of solution is from 0.5 to 5 ml. Accuracy attained is 80-90% by Lovibond tintometer 
and 70-80% by spectrophotometer. 

The reaction is specific for free sulphur and not influenced by compounds such 
as mercaptans, sulphides and saturated disulphides. Olefinic sulphur compounds, 
such as olefinic disulphides, undergo partial decomposition at the temperatures in- 


volved. 
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BOOK REVIEW 


THE CHEMISTRY AND TECHNOLOGY OF CEREALS AS FOOD AND FEED. Edited by 
Samuel A. Matz, Avi Publishing Co., Inc., Westport, Conn., 1959, 732 pp., $16. 


This book gives a survey of the chemistry and technology of the cereal industries. The book is directed 
to food technologists, agronomists, marketing and sales personnel, nutritionists and food chemists. 

The material is presented as a series of separate topics and each chapter is complete in itself. 
The book is divided into three sections. 

Section 1—The cereal grains: Wheat, corn, oats, barley, rye, sorghum, rice, millet, wild rice, 
adlay and rice grass. 

Section 2—Processing methods: Milling, commercial baking procedures, maccaroni production, 
manufacture of prepared mixes, wet milling, production of oils from cereal grains, feed manufacture, 
rice reese, manufacture of special dietary foods, malting, brewing, manufacture of breakfast 
cereals. 

Section 3—Functional and storage characteristics of cereal products: Characteristics of cereal 
starches, flavour staling resulting from lipid deterioration, adequacy of processed cereals in human 
nutrition, supplementation of cereal based animal feeds, effects and detection of insect and rodent 
infestation of cereals. 

The 25 chapters were written by 25 authors, all of them recognized authorities in their field. 
Though the editor stresses in the preface that breadth, rather than depth, of coverage was sought, 
the differences in coverage and presentation are too big to be explicable by the intention to cover 
so big a field in such small space. Whereas some chapters give a well presented and balanced sum- 
mary of the literature and of current processing methods, some chapters are actually a compilation 
of data, no effort being made to digest the material. The chapters of the first section are well handled 
and provide a convenient introduction into the botany, agriculture, handling and international 
trade of the main cereal grains. The major merit of the second section is in providing useful in- 
formation which is widely scattered over many technical and trade journals. Unfortunately, due 
to the large number of topics covered, some of the reviews are too thin to do justice to the subject 
being surveyed. With regard to the third section, there are given interesting reviews on cereal starches 
and on cereal lipids but for no obvious reason no chapters are devoted to proteins and to enzymes. 
The two chapters on nutrition give a fairly good coverage of the subject and provide a useful source 
of information. The treatment of insect and rodent infestation, though well written, is too thin to 
be of value to either the food technologist or laboratory worker. 


The book is well typed and illustrated, but some of the illustrations are not very descriptive and 
could be advantageously replaced by diagrams. 

In spite of these limitations, the book offers a wealth of material, provides valuable information, 
and the bibliography at the end of each chapter supplies ample information and references for the 
reader interested in a particular field. 
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